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محمدامین برفه‌ئی!! مهدا مرتضوی زنجانی 6910 عباس سروش" 
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چکیده در اين مقاله تأثیر جابه‌جایی گسل شیب‌لغز معکوس بر تونلی عمود بر خط گسل در پلان با نرمافزار اجزای محدود پلکسیس سه‌بعدی مطالعه 
شده و اثر زاویه گسل. نوع خاک منطقه و عمق قرارگیری تونل بر «مسیر انتشار گسلش». «نمودارهای جابه‌جایی سطح خاک» و «تغییر شکل و نیروهای 
وارد بر پوشش تونل» ارائه شده‌اند. نتایج تحلیلها نشان می‌دهد که تونل باعث ایجاد نواحی تمرکز کرنش در خاک زیر خود و تا حد زیادی مانع از رسیدن 
جابه‌جایی گسلش به لایه‌های خاک بالا سر خود (در محور تونل) شده حال آنکه موجب بالازدگی ذسبی خاک در دو طرف آن در سطح زمین می شود. 
علاوه بر آن در تونل عمیقتر نیروی وارده بر پوشش تونل (نیروی محوری و لنگر خمشی) در اثر جابه‌جایی گسل بیشتر است. همچنین با آنالیز حساسیت 
برای کاهش اثر مرز. نسبت ابعادی مناسب برای مدل‌سازی اندرکنش تونل و گسل پيشنهاد شده است. بررسیها نشان می‌دهد که نیروی محوری به خصوص 
در زوایای گسل ملایم به طول مدل حساس است حال آنکه لنگر حمشی در مدل سازی با ابعاد متعارف (1-۵ برابر عمق لایه حاک) به طور منا سب قابل 


ارزیابی می‌باشد. 


واژه‌های کلیدی انتشار گسلش. مدل‌سازی عددی اجزای محدود. تحلیل سه‌بعدی. توئل. پوشش بتنی. 
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(۱) کارشناسی ارشد. واحد دانشگاهی گرمسار, دانشگاه صنعتی امیرکبیر . 

(۲) نویسنده مسئول: استاد مدعو, واحد دانشگاهی گرمسار, دانشگاه صنعتی امیرکبیر. 2.6۴ ۵ ۵۷1۳ ۵۲)۵2 ۱۱۵۵۵۱۰۳۲ 
(۳) استاده دانشکده مهندسی عمران و محیط زیست. دانشگاه صنعتی امیرکبیر. 


۳۱ 


۳۲ 


مقدمه 

تا پیش از زمین‌لرزه‌های سال ۱۹۹۹ در ترکیه و تایوان (کوجایلی» 
دوزکه و چی‌چی) که موجب خرابیهای بسیار سازه‌های زیرزمینی 
گردید. باور عموم جامعه مهندسی بر این بود که سازه‌های 
زیرزمینی به دلیل محصور بودن در میان خاک از آسیبهای ناشی 
از لرزه در امانند. ولی بعد از این وقایع تحقیقات جدی و گسترده 
در مورد اثرات لرزش زمین و همچنین جابه‌جایی ماندگار ناشی 
از گسلش بر این سازه‌ها شدت گرفت. 

شید اتشار کسلش در لاه ماک عر سال ۱4۵۶ ترسط بر 
و همکاران [1] با روش اجزای محدود مورد بررسی قرار گرفت. 
آنها به تأثیرپذیری مسیر گسلش از زاویه صفحه گسل. انواع 
حرکت گسل. میزان جابه‌جایی گسل. عمق لایه خاک و نوع خاک 
محل پرداختند. در بررسیهای میدانی تأثیر نوع خاک بر گسلش؛ 
دریافتند که در خاکهای سختتر گسلش در جابه‌جایی پایه کمتری 
به سطح زمین می‌رسد. بعد از آن» بررسیهای متعدد آزمایشگاهی 
مانند یوهانسون و کاناگایی . آناستوپولوس و همکاران [2:3] و 
مدل‌سازیهای عددی مانند لین و همکاران. لاکیدیس و همکاران 
[4,5] تأثیر پارامترهای مختلف را بر پدیده گسلش در لایه خاک 
مورد مطالعه قرار دادند. در ادامه با مدل‌سازی عددی گسل 
معکوس در ماسه دولایه و چندلایه. وقوع انکسار در مسیر 
گسلش و شکل‌گیری الگوهای مختلف مسیر گسلش در خاک 
چند لایه نشان داده‌شد [6,7]. همچنین پدیده انتشار گسلش با 
روشهای عددی دیگر مانند روش المان مجزا [8-11] مدل‌سازی 
کل فرش وفتان دهاش شاک بت ار کسلین. اتلرمل] ۶ 
تأثیر عمق دفن در اندرکنش گسل و فونداسیون [14] ابعاد بیشتری 
از پدیده گسلش را روشن کرد. 

در سالهای اخیر بررسی اندرکنش گسل و سازه‌های مدفون 
مانند لوله و تونل با روشهای آزمایشگاهی انجام شده است. 
بازیار و همکاران [15] با بررسی اندرکنش گسل شیب لغز 
معکوس و تونل با کمک آزمایش سانتریفیوژ دریافتند سازه‌های 
زیرزمینی نظیر تونلها می‌توانند بر مسیر گسلش تأثیر بگذارند که 
با افزايش صلبیت تونل و افزايش عمق قرارگیری آن فشار ناشی 
از گسلش بر تونل بیشتر می‌گردد. در پژوهش آنها راستای طولی 
تونل موازی خط گسل در پلان قرار داشت. 

رجحانی و همکاران (۲۰۱۲) [16] با بررسی نتایج چهار 
آزمایش سانتریفیوژ با تمرکز بر رفتار لوله‌های فولادی پیوسته 


نشریه مهندسی عمران فردوسی 


مدلسازی عددی سه بعدی آند رکن شگسل شیب لغز و توتل... 


مدفون تحت اثر گسلش معکوس. نشان دادند با تغییرات قطر و 
عمق دفن لوله‌هاه مکانیزم تغییر شکل لوله و نوع خرابی به طور 
اساسی تغییر می‌کند. کیانی و همکاران [17] با انجام آزمایشهای 
سانتریفیوژ, تأثیر انتشار گسلش نرمال را بر روی تونل با پوشش 
قطعدا تسین گنل 

پررسی پدیده گسلش در اندرکنش با تونلها با روشهای 
عددی هم صورت گرفته است. جوشی و همکاران [18] با روش 
اجزای محدود و فرض ساده المان تير برای لوله. تأثیر زاویه 
قرارگیری خط لوله نسبت به خحط گسل را بررسی کردند و 
همچنین به راهکارهای کاهش خطر گسلش همچون کاهش عمق 
دفن در محل گسل و استفاده از خاکریزی با مصالح سست اشاره 
کردند. هاشم‌زاده و سروش [19 با مدل‌سازی سه‌بعدی لوله 
مدفون تحت گسل امتدادلغز نشان دادند کرنش طولی خط لوله 
در ماسه‌های سختتر افزایش می‌یابد و کمانش موضعی لوله در 
جابه‌جایی گسل کمتری به وقوع می‌پیوندد. 

راششریا فاشکاران (20] ۵ تسش فا یرو یرای 
تونلهای بتنی کم‌عمق را تحت اثر گسلش شیب لغز معکوس و 
تما فطالمه کر ور شام تداوزن کهبا آفر ای ناسا 
پیشتری در هر دو نوع گسل به تونل وارد می‌شود. 

در مقاله حاضر به کمک مدل‌سازی عددی با نرم‌افزار اجزای 
محدود پلکسیس سه‌بعدی [21] به تحلیل سه‌بعدی گسل شیب لغز 
معکوس با راستای عمود بر محور طولی تونل پرداخته شده و در 
آن اثر زاویه گسل, نوع خاک و عمق قرارگیری تونل بر مسیر 
انتشار گسلش, جابه‌جایی سطح زمین. تغییر شکل پوشش تونل 
و نیروهای ایجاد شده در پوشش تونل بررسی شده است. لازم به 
ذکر است که منظور از نیروی محوری و لنگر خمشی در اینجا؛ 
نیروهای ایجاد شده در پوشش تونل (8فتل) می‌باشد. همچنین 
در این مقاله تلاش شده است یافته‌های عددی به نتایج کاربردی 
منجر شود و توصیه‌هایی برای طراحی ایمن تونل در برابر گسلش 


ان یگ 


معرفی کلی مسئله 
برای مدل‌سازی عددی سه بعدی انتشار گسلش شیب لغز 
معکوس» حجمی از خاک به عمق ۲۰ متر و به طول چهار برابر 
عمق لایه خاک (۸۰ متر)؛ به منظور کاهش اثر مرزهاه و نیز عرض 


سال سی و ششم شماره یک, ۱۶۰۲ 


محمدامین برفه‌شی - مهدا مرنضوی زنجانی - عباس سروش 


۰ متر در نظر گرفته شد که به جهت تقارن عرضی مسئله نیمی 
از عرض به اندازه ۲۰ متر مدل شده است. سطح خاک افقی در 
نظر گرفته شده و حرکت گسل که از اعماق سنگ بستر صلب 
شروع شده با زاویه 0 به فصل مشترک سنگ بستر خاک می‌رسد. 
برای مدل‌سازی حرکت گسل» قسمت سمت راست مدل 
(فرودیواره) ثابت و جابه‌جایی به سطح پایین سمت چپ و 
قسمت عمودی سمت چپ (فرادیواره) مدل وارد می‌گردد (شکل 
فاصله گسل از سطح قائم سمت فرادیواره ۳۰ متر و زاویه 
گسل با افق به دو صورت ۷۰ درجه و قائم در نظر گرفته شده 
است. میزان جابه‌ جایی سنگ بستر ۰,۸ متر در راستای قائم فرض 
شده است که مولفه افقی جابه‌جایی با توجه به زاویه گسل 
مشخص خواهد شد. جابه‌جایی به صورت استاتیکی بر بدنه خاک 
اعمال می‌شود و از آثار دینامیکی حرکت زمین صرف نظر 
می‌گردد. 

تونلی بتنی به شعاع ۲/۱۲ متر در دو موقعیت سطحیتر و 
عمیقتر (عمق قرارگیری محور تونل از سسطح خاک به ترتیب ٩‏ 


متر و ۱۱ متر) در نظر گرفته شده است. 
۱ 


شکل ۱ هندسه سه بعدی مدل‌سازی اجزای محدود 


لایه خاک در نظر گرفته شده در این مسئله از نوع ماسه 
شک بوده و دو نوع ما سه متراکم (5) و سست (,1) برای 
مدل‌سازی در نظر گرفته شده است. مشخصات دو نوع خاک در 
جدول (۱) آمده است. مشخصات تونل و مصالح آن نیز به شرح 
جدول (۲) می‌باشد. 
بدین ترتیب حالتهای مختلف در نظر گرفته شده در این تحقیق 


با تغییر موارد زیر مشتمل بر ۱۲ حالت می‌باشد: نوع خاک (ماسه 


سال سی و ششم شمارة یک, ۱۶۰۲ 


۳۳ 


متراکم و سست) (جدول ۱ زاویه گسل (1۰ درجه و قائم) 
شرایط میدان آزاد (عدم حضور تونل) و قرارگیری تونل در دو 
موقعیت سطحیتر و عمیقتر (عمق قرارگیری محور تونل از سطح 
خاک به ترتیب ۱ متر و ۱۱ متر) (جدول ۲) 


حدول ۱ مشخصات خاک 


جدول ۲ مشخصات بتن پوشش تونل 


رززن ۳۴ ۳۳/۵ ۲۵ ۲ | الاستیک 


مدل‌سازی عددی. مدل موهر-کولمب مدل مناسبی برای بررسی 
مدل استفاده شده است. در اینجا برای معرفی رفتار خاک مدل 
رفتاری موهر-کولمب الاستیک-کاملاً پلاستیک با قانون جریان 
غیرهمبسته (1۲016 110۷ 002890012]60) در نظر گرفته شد. 
همچنین در مورد پوشش بتنی توئل فرض شده است که رفتار 
الاستیک داشته باشد. 

برای مدل‌سازی فصل مشترک خاک و بتن از المانهای رابط 
(101611200) دارای اصطکاک و چسبندگی استفاده می‌شود. 
روابط زير برای محاسبه مقدار زاویه اصطکاک و چسبندگی آنها 


اواته کف آبگ: 


)۱( م ق مررک ‏ م م2 20 
۳( سجن ویک روز 
با توجه به رابطه بالا و ,2 پٍ شنهادی بازیار و همکاران 
ر دعر هلان 7 1داعدون به ناست هی ای 

نوع المان مورد استفاده برای خاک المان چهاروجهی ده 
گرهی و المان مورد استفاده برای پوشش تونل المان 01216 مثلثئی 


شتی کرهی استتة 


نشريهة مهندسی عمران فردوسی 


۳۶ 


مدلسازی عددی سه بعدی اند رکن شگسل شیب لغز و توتل... 


7۱/۹ 
۱۵۵ 
۵۲۵۲۵۵۲۵ 
۱ 


شکل ۲ نمایش مش‌بندی و تغییر شکل زمین و تونل بعد از وقوع گسلش 


به منظور صحت‌سنجی مدل عددی آزمایش سانتریفیوژ 
بازیار و همکاران [15] به صورت عددی مدل‌سازی شد و 
جابه‌جایی گسل شیب‌لغز معکوس با زاویه گسل ٩۰‏ درجه با سه 
اندازه جابه‌جایی قائم مختلف سنگ بستر در شرایط بدون حضور 
توئل و با حضور تونل اعمال شد. مقایسه نتایج تغییر شکل سطح 
زمین در آزمایش سانتریفیوژ با مدل‌سازی عددی انطباق خوبی 
را نشان می‌دهد. شرایط مش‌بندی و تغییر شکل مسئله در حالت 
گسل با زاویه قائم و تونل سطحی‌تر در شکل (۲) نشان داده شده 


است. 


بر رسی نتایج 


سیر انتشا رگسلش (جابه‌جایی زمین و تونل). عملکرد گسلها 
عموماً موجب به‌وجود آمدن جابه‌جایی‌های برشی مهمی در 
سطح زمین و همچنین سازه‌های سطحی و مدفون می‌شود. این 
جابه‌جاییهای اضافی تحمیل شده بر سازه, عملاً در طراحیهای 
اولیه در نظر گرفته نمی‌شود و ممکن است پایداری سازه‌ها را 
برهم بزند؛ لذا بررسی این آثار بر روی سازه‌ها بسیار با اهمیت 
تشن 

شور ان پزوهتن سیر انقشار کسلش: که کرنتن برشی مه 
معرف آن است در حالتهای مختلف مسئله و برای دو زاویه 
مختلف گسل (زاویه 1۰ و ٩۰‏ درجه)؛ دو نوع خاک (خاک ماسه- 
ای متراکم و یا سست) و تونل سطحیتر و عمیقتر بررسی شده و 
به نحوه اثر تونل بر این یدیده پرداخته تال م: ات از شکل (۳) 
نمونه‌ای اوقفی او سای را در دو حالت تونل سطحیتر 


سست و گسل با زاویه ۰ درجه در مدل سه بعدی نشان می‌دهد. 


نشريه مهندسی عمران فردوسی 


در این شکل وضعیت توزیع کرنش در مقاطع مختلف نشان داده 
شده است. 

شک (۳ال) برضی سبه بعتی از وسط توت وعاک است: 
صفحه قائم میانی که تونل در آن قرار دارد شکل (۳سب) و 
ی ی ی 
تضنویر عناسبی, از تشیرشگلها ازانه می‌دهته. ور این شکل 
مقششصی. اسک که تسیر گسلکن کن مسلوهه کون با کبارن‌ها 
متفاوت است (شکل "سب در مقایسه با شکل ۳-ج) . این 
تفاوت در رخنمون گسلش در سطح زمین در بالای محور تونل 
و کناره‌ها هم مشهود است (شکل ۳۲-د). 

بررسی این شکل نشان می‌دهد که تونل باعث ایجاد نواحی 
تمرکز کرنش در خاک زیر خود و تا حد زیادی مانع از رسیدن 
جابه‌جایی گسلش به لایه‌های خاک بالاسر خود شده است 
(شکل ۳الف و ب). در حالی که با فاصله گرفتن از تونل در 
جهت عرضی (در جهت محور 6 تأثیر حضور تونل نیز تقریبً 
از بین رفته است (شکل ۳-د) و در صفحات کناری (شکل ۳- 
ج) مسیر گسلش مطابق شکل زمین آزاد (15610 8506) شکل 
میگیرد. 

به بیان ساده‌تر حضور تونل باعث ایجاد بالازدگی نسبی در 
خاک دو طرف می‌شود. توجه به این تغییر شکل غیریکنواخت؛ 
به لحاظ مهند سی اهمیت دارد زیرا سازه‌های سطحی را تحت 
تغییر شکل غیریکنواخت قرار می‌دهد و به خصوص در بالا سر 
سازه‌های مدفون دارای صلبیت زیاد (مانند تونلهای شهری یا 
لوله‌های قطور) خطرساز خواهد شد. 


سال سی و ششم» شمارهة یک, ۱2۰۲ 


محمدامین برفه‌ثی - مهدا مرنضوی زنجانی - عباس سروش 
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۳۵ 


بالازدگی خاک در کناره‌ها 


۷ 
تمرکز کرنش زیر توفل 


نرسیدن گسلش به 


سطح در بالاسر نز 0 


شکل ۳ بررسی نحوه انتشار گسلش (کانتورهای کرنش برشی بیشینه) در دو حالت تونل سطحیتر (سمت راست) و توئل عمیقتر (سمت چپ). با فرض خاک 


مسست و گسل با زاویه 1۰ درجه: الف) تصوير سه بعدی» ب) برش از محور تونل صفحه 20 ج) تصوير در صفحه کناری. ۷20 د) تصویر در سطح زمین 


اثر سه عامل مهم یعنی جنس خاک زاویه گسل و عمق 
زار گنر کرت در اکن بر کی شکلای تاشی از العقای سل 
به صورت نواحی کرنش برشی بیشینه. در شکل (4) نشان داده 
شده است. مقایسه کرنشها در تونل عمیق و کم‌عمق برای خاک 
یکسان (مثلاً شکل ۶-الف و پ برای خاک متراکم) نشان می‌دهد 
که هنگامی که تونل عمیقتر باشد. بیشتر جابه‌جایی ناشی از 
گسلش را در خود جذب می‌کند و تمرکز کرنش کمتری در سطح 
زمین بالاسر تونل و کناره‌ها مشاهده می‌شود. از سوی دیگر 
هنگامی که عمق تونل بیشتر باشد. در زیر تونل و پوشش آن 
کرنشهای بیشتری توسعه می‌یابد که منجر به تنشهای بیشتر 
خواهد شد. این مطلب در ادامه در بررسی نیروی محوری و لنگر 
خمشی در تونل کم عمق و عمیق (شکل ۱۱ و ۱۲) خواهد آمد. 

در خاک متراکم نسبت به خاک سست با عمق تونل یکسان 
(مثلاً مقایسه شکل 4-لف و ث) تمرکز کرنش بیشتری در زیر 
توتل و سطح خاک مشاهده می‌شسود. این مطلب با یافته‌های 
پیشین در مورد تمرکز بیشتر کرنش در خاکهای متراکم. که به 
معنی محدوده باریکتر گسلش است. سازگار است [22]. 

تحقیقات پیشین انتشار گسلش در بررسی جابه‌جاییهای 
سطح زمین در شرایط لایه خاک بدون حضور تونل ۳۲۵6) 
(۳1610 نشان داده است که افزایش تراکم خاک و کمتر شدن 


سال سی و ششم, شماره یک, ۱۶۰۲ 


زاو به تن (0) هر دو باعت انتقال جابه‌جاییها به سسمت 
فرودیواره ی گر کل (مرتضوی و همکاران. 0۳۰۲ [23]. این 
مطلب در پژوهش حاضر هم مورد بررسی قرار گرفت و نتایج 
(شکل ۵) یافته‌های قبلی را تأیید کرد. 

شکل (1) تغییر شکل سطح زمین را در خاک متراکم و زاویه 
گنسلی ۰ درجه در دو حالت تونل عمیقتر و تونل سطحیتر در 
مقایسه با حالت زمین آزاد (عدم حضور تونل) نشان می‌دهد. 
مشاهده می‌شود که حضور تونل باعث می‌شود تغییر شکل سطح 
زمین در بالاسر آن توزیع ملایمتری (شیب ملایمتر) داشته باشد 
که در توئل عیق: این تایر تین پیشعر می‌شود. در واقع حضور 
توئل و هر چه عمیقتر بودن آن» موجب تمرکز کرنش در زیر خود 
و کاهش تمرکز کرنش در لایه‌های بالاسر خاک (محور تونل) 
می‌شود. جنان چه در این شکل مشاهده می‌شود در سطح زمین 
جابه‌جایی ۰/۹ متر از گسل؛ در شرایط عدم حضور گس در 
عرض ۱۱ متر و در شرایط حضور تونل عمیق در ۱٩‏ متر پخش 
شده است. این امر در ترکیب با بالازدگی خاک در دو طرف تونل 
(شکل ۳) وضعیت غیریکنواختی از نظر توزیع عرضی تغییرشکلها 
در سطح زمین ایجاد می‌کند. 


نشریه مهندسی عمران فردوسی 


۳۹ مدلسازی عددی سه بعدی اند رکن شگسل شیب‌لغز و تونل... 


شکل ۶ کانتورهای کرنش برشی بیشینه در اثر گسلش در لایه خاک: الف) خاک متراکم. زاویه ۰ درجه و تونل سطحیتر ب) خاک متراکم زاویه ۰ درجه و 
تونل سطحیت پ) خاک متراکم زاویه ۰ درجه و تونل عمیقتر ت) خاک متراکم زاویه ۰ درجه و تونل عمیقتس ث) خاک سست. زاویه 7۰ درجه و تونل 


سطحیتر. ج) خاک سست. زاویه ۰ درجه و تونل سطحیتر» چ) خاک سست. زاویه ۰ درجه و تونل عمیقت ح) خاک سست. زاویه ۰ درجه و تونل عمیقتر 


خاک متراکم. زاویه گل ۶۰ سس 


90 80 7 60 50 40 30 20 10 0 
طول نمونه (17) 


شکل ۵ نمودار تغییر شکل سطح زمین بعد از گسلش در حالت میدان آزاد (عدم حضور تونل) با تغییر نوع خاک و زاویه گسل 


نشریه مهندسی عمران فردوسی سال سی و ششم شمار یک. ۱۶۰۲ 


محمدامین برفه‌ثی - مها مرتضوی زنجانی - عباس سروش ۳۷ 


[ م] 


توئل عمیق سس 
ترئل سطحی سس 
عدم حضور توئل- - - 


90 80 70 60 50 40 30 20 10 0 
طول نمونه (0) 
شکل ٩‏ تغییر شکل سطح زمین در خاک متراکم و زاویه گسل ٩۰‏ درجه (توزیع ملایمتر شیب زمین در حالت بدون تونل) 


[ یا 


ك 

۸ 

4 

ك 

خاک متراکم با زاویه گسل ۶۰ درجه سس : 
خاک متراکم با زاوبه گسل قائم ۳ 1 
خاک سست با زاوبه گسل ۶۰ درچه- << 1 
خاک سست با زاویه گسل قائم -۰- ِ 


90 


شکل ۷ تأثیر جنس خاک بر میزان جابه‌جایی قائم ایجاد شده در سقف پوشش تونل (تونل عمیقتر) در دو زاویه قائم و ۱۰ درجه 


اس 


8 

۸ 

ب ۳ 

۱ ضِ 

خاک متراکم. تونل عمیق سس 3 

خاک متراکم. تونل سطحی وه ۴ 
خاک سست. تونل عمیق یت 4 

مت مکی سب 


90 


طول نمونه (00) 


شکل ۸ تأثیر عمق تونل بر جابه‌جایی قائم سقف پوشش تونل در زاویه گسل ۰" درجه در دو نوع خاک متراکم و سست 


سال سی و ششم شماره یک, ۱۶۰۲ نشریه مهندسی عمران فردوسی 


۳۸ 


شکل (۷) میزان جابه‌جایی قائم ایجاد شده در سقف پوشش 
تونل بر اثر گسلش در دو نوع خاک متراکم و سست و دو زاویه 
قائم و 1۰ درجه در حالت حضور تونل عمیقتر را نشان می‌دهد. 
شکل (۸) جابه‌جایی قائم سقف پوشش تونل در دو نوع خاک 
متراکم و سست در زاویه گسل ٩۰‏ درجه را برای عمقهای مختلف 
تونل نشان می‌دهد. از این دو شکل می‌توان نتیجه گرفت که تراکم 
بیشتر خاک و ملایمتر شدن زاویه صفحه گسل هر دوء باعث 
انتقال جابه‌جایی‌های ایجاد شده در پوشش تونل به سمت 
فرودیواره می‌شوند. شیب ملایمتر جابه‌جایی در شرایط تونل 


عمیق در این شکل هم دیده می‌شود. 


نیروهای وارد بر پوشش تونل. در اثر جابه‌جایی گسل. پوشش 
تونل تحت نیرو و لنگر در جهتهای مختلف قرار می‌گیرد. در 
نرم‌افزار پلکسیس سهبعدی, داده‌های خروجی برای یک صفحه 
فان شایل عفی‌شکلها و تروها مطانق شکل (4) می‌باشا 
پدیده گسلش به طور عمده باعث به وجود آمدن خمش در 
راستای طولی تونل (۷111) و نیروهای کششی (فشاری) ۱1 در 
طول محور تونل می‌شود. طبق قرارداد در پلکسیس نیروی 
محوری در کشش مثبت (+) و در فشار منفی (-) در نظر گرفته 
می‌شود. علامت ممانهای خمشی و نیروهای برشی نیز به سیستم 
مختصات محلی صفحه بستگی دارد. 

بررسی نیروی محوری ایجاد شده در ناحیه سقف پوشش 
تونل در طول محور آن در گسلهای با زاویه ۶0 ۸۰ ۷۵ و ٩۰‏ 
درجه برای خاک متراکم و تونل سطحی (شکل ٩۰‏ نشان می‌دهد 
که با ملایمتر شدن زاویه گسل» نیروی محوری فشاری ایجاد شده 
در پوشش تونل به طور چشمگیری افزايش می‌يابد. در اینجا 
زوایای گسل 1۵ و ۷۵ درجه برای بررسی کاملتر به تحلیلها اضافه 
شده‌است. 

علاوه بر اين. شکل (۱۰) نشان می‌دهد که در گسل قائم 
بخشی از سقف تونل به کشش می‌افتد ولی در گسلهای با زاویه 
0 تا ۷۵ درجه سقف تونل تماما تحت فشار قرار می‌گیرد. 

همچنین از بررسی تأثیر عمق تونل در زاویه‌های گسل قائم 
و ۰ درجه در حالت خاک متراکم (شکل ۱ به نظر می‌رسد که 


نشریه مهندسی عمران فردوسی 


مدلسازی عددی سه بعدی اند رکن شگسل شیب‌لغز و تونل... 


در تونل سطحیت محدوده تغییرات نیروی محوری کمتر است 
(به عنوان مثال در زاویه ٩۰‏ درجه. تفاضل ماکزیمم تا مینیمم 
نیروی محوری در تونل سطحیتر ۲۱۰۰ کیلونیوتون و در تونل 
عمیقتر ۳۰۰۰ کیلونیوتون است). این بدان معنی است که در 
تنشهای محدود کننده بیشتر (که در عمق بیشتر اتفاق می‌افتد) 
نیروی محوری بیشتری در اثر گسلش در تونل شکل می‌گیرد که 
در طراحی باید مورد توجه قرار گیرد. 

شکل (۱۲) نمودار لنگر خمشی را در حالت خاک متراکم 
در دو زاویه 1۰ درجه و ٩۰‏ درجه و در دو حالت تونل سطحیتر 
و تونل عمیقتر نشان می‌دهد. براساس این نتایج در گسل با 
زاویه ٩۰‏ درجه ذسبت به گسل قائی نمودار لنگر خم‌شی ایجاد 
شده در سقف توئل به سمت فرودیواره جابه‌جا شده‌است و در 
تونل عمیقتر مقادیر لنگر خمشی بیشتر از تونل سسطحیتر (در 
زاویه گسل عشابه) می‌باشد. این مسئله بار دیگر اهمیت طراخی 


پوشش تونل در تونلهای عمبة را نشان می‌دهد. 


شکل ٩‏ معرفی نیروها و لنگرها در پوشش تونل 


سال سی و ششم, شماره یک, ۱۶۰۳ 


محمدامین برفه‌شی - مهدا مرنضوی زنجابی - عباس سروش 


یا 


ح‌‌ 


گسل با زاوبه ۴۵ درچه سس 
گسل با زاویه ۶۰ درچه س.... 
گسل با زاویه ۷۵ درجه تسد 
گسل با زاویه ٩۰‏ درجه 7 


ی و دا 
1 


۳۹ 


هه -عحی ی 10000 
2 »۰ 2000 


-30000 
-40000 


نیروی محوری (10) ۲1 


-000 
-0000 
-0000 
-0000 


شکل ۷۰ نیروی محوری سقف تونل در طول محور تونئل در گسلهای با زاویه ۵ ۵ و ۰ درجه در خاک متراکم و تونل سطحی 


ونل عیق‌تره زاویه گسل 60 سب 


تونل عمی‌ترهزاویه گسل 90 یچ 
تونل سطحی‌ترهزاوبه گم 60 - -- و ۳ 
تونل سطحیتره زاویه گسل 90 - - - # سیم 
۳ 2و6 4 ۷ 2 
کم ۱ لك 
ی ۳ 


7 00۷ 


صول نموئه (۱0) 


شکل ۱۱ تأثیر عمق تونل بر نیروی محوری در طول مدل در زاویه‌های 
گسل قائم و ۰ درجه در خاک متراکم 


تونل ععیقتر, زاویه گسل ۶۰ درچه سس 

تونل ععیقتر, زاویه گسل ٩۰‏ درجه 
توئل سطحی‌تر, زاوبه گسل ۶۰ درجه - - - ۳ 40 
تونل سطحیتر, زاویه گسل ٩۰‏ درجه .-- 


نگر خمشی (۱۲.۱) 3/11 


طول نمونه (۳) 
شکل ۱۲ : نمودار لنگر خمشی در حالت خاک متراکم در دو زاویه 7۰ در 
جه و ٩۰‏ درجه و در دو حالت تونل سطحیتر و تونل عمیقتر 


سال سی و ششم شمارة یک, ۱۶۰۲ 


بررسی اثر مرز 

در مدل‌سازی سازه‌های طویل همچون تونل, افزایش طول مدل 
به افزايش تلاش محا سباتی منجر می شود. لذا لازم است طول 
بهینه برای مدل‌سازی سازه انتخاب شود. معیار اين انتخاب؛ 
حصول اطمینان از صحت نتایج است. 

در مدل‌سازی اولیه با توجه به توصیه مراجع برای لایه خاک 
طول مدل ۸۰ متر (؛ برابر عمق لایه خاک) انتخاب شد. به منظور 
بررسی اثرات مرزبندی, با افزايش طول مدل از ۸۰ متر به دو 
طول ۱۸۰ متر (اضافه کردن ۵٩۰‏ متر به ابتدا و انتهای مدل) و ۲۸۰ 
متر (اضافه کردن ۱۰۰ متر به ابتدا و انتهای مدل)» مورد بررسی 
قرار گرفته است. شایان ذکر است که به عنوان معیار بررسی از 
بین نتایج متنوع می‌توان پارامترهای هدف متعددی را؛ مانند 
تغییر شکل سطح زمین. جابه‌جایی و کرنش در سقف و کف 
تونل» و نیروها در پوشش تونل. مد نظر قرار داد. از این بین؛ 
تنشها (نیروها) در پوشش تونل می‌توانند نماینده مناسبی برای 
سایر پارامترها باشند زیرا به دلیل رفتار الاستیک بتن, تنشها و 
کرنشها به صورت خطی با یکدیگر تناسب دارند. از سوی دیگر 
بررسی مقادیر نیروها و لنگرها. امکان انجام طراحی سازه‌ای را 
هم فراهم می‌کند. لذا در اين مقاله برای مقایسه اثر مرز در 
مدلهای با طول ۰۸۰ ۱۸۰۰۲۸۰ از نمودارهای نیروی محوری 
(11) و لنگر خمشی (۷11) استفاده شده است: 

محل اعمال سل در تمامی مدلها بر روی مبداً مخد صات 
قرار دارد. در شکل (۱۳) نمودار نیروی محوری (۷1) وارد بر 
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پوشش توئل در طول مدل. در مدلهای با طولهای ۸۰ ۱۸۰ و 
۰ متری و برای زاویه‌های گسل ۰4۵ 1۰ و ٩۰‏ درجه رسم 
شده است. در این شکل نتایج مدلهای ۸۰ ۱۸۰ و ۲۸۰ متری. 
از تفاوت در طول آنها قابل تشخیص می‌باشند. 

با توجه به شکل (۱۳) در زاویه گسل کمتر (گسل 1۵ 
درجه در مقایسه با گسل ٩۰‏ درجه» حساسیت نیروی محوری 
به طول مدل انتخابی زیاد است. به طوری که برای زاویه سل 
0 درجه. با افزايش طول مدل از ۸۰ متر به ۲۸۰ متر (واقعیتر 
شدن مدل سازی)» نیروی محوری فشاری وارد بر پوشش تونل 
کاهکن قایل توجهی می‌باید. این در حالی است که قر سل با 
زاویه قائم با افزايش طول مدل. تغییر چندانی در نیروی فشاری 
وارد بر پوشش تونل (1) ایجاد نمی‌شود. 

بنابراین در مدل‌سازی انتشار گساش در زوایای گسل 
ملایمتر. جهعت حذف اثرات نامطلوب مرز, انتخاب طول 
بیشتری برای مدل توصیه می‌شود. چنان چه در مدل‌سازی 
عددی. ابعاد متداول در ادبیات (طول حدود ۵-۶ برابر عمق 
لایه خاک) | ستفاده شده با شد. تسیر نتایج نمودارهای نیروی 


محوری, نیاز به اصلاح در زوایای گسل کمتر (به خصوص 1۵ 
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نمونه‌های با زاویه گسل ٩۰‏ درجه تس 
تمونه‌های با زاوبه گسل ۶۰ درچه سب 
تمونه‌های با زاوبه گسل ۴۵ درجه -.- 


مدلسازی عددی سه بعدی اند رکن شگسل شیب‌لغز و تونل... 


درجه و ٩۰‏ درجه) و کاهش مقدار نیروی محوری دارد (مانند 
شکلهای ۱۱ و ۱۲ اين مقاله). 

در شکل (۱۶) نمودار لنگر خمشی (۷۲11) وارد بر پوشش 
تونئل در طول مدل. در مدل با طولهای ۸۰ ۱۸۰ و ۲۸۰ متری و 
برای زاویه‌های گسل 0 درجه رسم شده‌اند. همان طور که 
م شاهده می شود. افزایش طول مدل تغییر مح سو سی در لنگر 
خمشی وارد بر پوشش تونل ایجاد نمی‌کند. 

بررسیها نشان می‌دهد که در زوایای گسل کمتر که نسبت 
مولفه افقی جابه‌جایی به ملفه قائم آن بیشتر است. توجه به اثر 
مرزها در مدل سازی (فا صله مرزها از محل گ سلش يا ذ سبت 
طول به عمق مدل) برای برآورد صحیح نیروی افقی وارد بر 
پوشش تونل. و در نتیجه تغییرشکلها و کرنشهای محوری؛ 
ضرورت دارد. 

بدین ترتیب در تخمین نیروی محوری واقعی وارد بر 
پوشش تونل در مسئله انتشار گسلش برای زوایای گسل کمتر 
جهت حلذف اثرات نامطلوب مرز. مدل طویلتر ( سبت طول به 
عمق ۱۵-۱۰). 
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20000 
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-60000 


-0000 


طول نموئه (0) 


شکل ۱۳ نمودار نیروی محوری وارد بر پوشش تونل (1) در طول مدل. در مدلهای با طولهای ۰ ۰ و ۲۸۰ متر و برای سه زاویه ۶۵ ٩۰‏ و ٩۰‏ درجه 
در شرایط خاک متراکم و تونل عمیقتر (مقادیر مثبت کشش و منفی فشار) 
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شکل ۱۶ نمودار لنگر خمشی وارد بر پوشش تونل (۷]11) در طول مدل. در مدلهای با زاویه گسل 4۵ درجه و با طولهای ۸۰ ۱۸۰ و ۲۸۰ متر 


در شرایط خاک متراکم و تونل عمیقتر 


پيشنهاد می‌شود در حالی که برای زوایای نزدیک به قائم 
نسبت طول به عمق ۶ کفایت می‌کند. همچنین برای تخمین لنگر 
خمشی وارد بر پوشش تونل, و تغییرشکلها و کرنشهای خمشی؛ 
در تمام زوایای گسل نسبت طول به عمق ۵-6 مناسب می‌باشد. 


طراحی سازه‌ای پوشش بتن در مقابل گسلش 


بررسی عملکرد سازه بتنی در برابر نیروهای وارده برمبنای 
آیین‌نامه بتن ایران (مبحث نهم مقررات ملی ساختمان) [24]. 
برای یک مقطع از سازه بتنی با طول واحد و با نیروهای خروجی 
مدل‌سازی عددی انجام می‌شود. ابعاد مقطع یک مستطیل به طول 
واحد یعنی ۱۰۰۰ میلیمتر و عرض ۲۶۰ میلیمتر که ضخامت بتن 
سازه بتنی می‌باشد. در نظر گرفته شده و فاصله آرماتورها مطابق 
استانداردهای اجرایی برابر ۷۵ میلیمتر در نظر گرفته می‌شود. در 
شکل (۱۵) هندسه مقطع رسم شتلام: اسجنت: 

در اینجا پوشش تونل تحت اثر توام نیروی محوری و خمش 
قرار گرفته است. لذا لازم است اندرکنش نیروی محوری و خمش 
مطابق روشهای مرسوم طراحی برای آن در نظر گرفته شود. اما 
به دلیلی که در ادامه توضیح داده خواهد شد. ابتدا ماکزیمم 
ظرفیت خمشی خالص (بدون حضور نیروی محوری) و ماکزیمم 
نیروی محوری (بدون حضور لنگر خمشی) مقطع فوق را به 
دست می‌آوریم. 
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شکل ۱۵ هندسه عرض واحد از پوشش بتنی تونل 
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۳۲ 


بنابراین حداکثر ظرفیت خمشی پوشسش تونل ۱۷۹ 
کیلونیوتون-متر می‌باشد که از ماکزیمم مقادیر لنگر خمشی وارد 
بر پوشش تونل (۳۰ کیلونیوتون مطابق شکل ۱۶) بیشتر است. 

در ادامه حداکثر ظرفیت فشاری و کششی مقطع بتنی را به 
فک ی انا 


تحمل کششی این مقطع (دو ردیف آرماتور بالا و پایین): 


)۵ ۷ 2 ور(2)090 < 7 


تحمل فشار این مقطع: 


6 < 06 << 7 0۰( 


همان طور که دیده می‌شود با توجه به مقادیر نیروی محوری 
در شکل (۱۳) (ماکزیمم فشار و کشش ۲۰۰۰۰ کیلونیوتون در 
گسل ٩۰‏ درجه و ماکزیمم فشار ۳۸۰۰۰ کیلونیوتون در گسل 10 
درجه). مقطع بتنی فوق. علی رغم ظرفیت تحمل خمش قادر به 
تحمل چنین نیروی محوری (حتی بدون حضور لنگر خمشی) 
نمی‌باشد. از این رو راهکارهای جایگزین کاهش خطر گسلش 
مانند استفاده از اتصالات منعطف پیشنهاد می‌گردد [124 


جمع‌بندی و نتیجه گیری 
در این مقاله اثر پارامترهای مختلف شامل نوع خاک زاویه گسل 
و عمق تونل بر مسیر انتشار گسلش. نمودارهای جابه‌جایی سطح 
خاک و همچنین تغییر شکل و نیروهای وارد بر پوشش تونل 
بررسی شد و ظرفیت خمشی و محوری پوشش تونل برای تحمل 
مقادیر حاصل بررسی شد. این نتایج را می‌توان به صورت زیر 
بیان کرد: 
الف) تأثیر تونل در تغییر شکل خاک: 
تونل باعث ایجاد نواحی تمرکز کرنش در خاک زیر خود و 
تا حد زیادی مانع از رسیدن جابه‌جایی گسلش به لایه‌های 
خاک بالاسر خود می‌شود. 
حضور تونل باعث ایجاد بالازدگی نسبی در خاک دو طرف 


می‌شود. توجه به تفاوت تغییر شکل سطح زمین در بالاسر 
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تونل و کناره‌هاه به احاظ مهندسی اهمیت دارد زیرا 
سازه‌های سطحی را تحت تغییر شکل غیر یکنواخت قرار 
می‌دهد. 

حضور تونل موجب توزیع ملایمتر جابه‌جاییها در سطح 
زمین (کاهش گرادیان حداکثر تغییر شکل سطح زمین) در 
بالاسر خود می‌شود که با افزايش عمق تونل این تأثیرگذاری 
ناحیه برشی تمرکز يافته در زیر تونل به سمت فرودیواره 
توسعه می‌یابد. لذا نسبت به شرایط زمین آزاد. مناطق 


فرودیواره مورد تهدید بیشتری هستند. 


ب) نیروهای ایجاد شده در پوشش تونل نانز کسبات ره 


در گسلهای با زاویه ملایمتر نیروی محوری وارد بر پوشش 
تونل بیشتر می‌باشد. در گسل نزدیک به قائم بخشی از سقف 
تونل به کشش می‌افتد ولی در گسلهای با زاویه کمتر سقف 
تونل تماما تحت فشار قرار می‌گیرد. 

در تونل عمیقتر نسبت به تونل سطحیتر» نیروی محوری و 
لنگر خمشی وارد بر پوشش تونل افزايش می‌یابد. این 
مطلب اهمیت زیادی در طراحی پوشش تونلهای عمیقتر 
دارد که باید مورد توجه طراحان قرار گیرد. 

در مدل‌سازی اندرکنش گسل و تونل در زوایای گسل 
نزدیک به قائم. انتخاب طول ٩-۶‏ برابر عمق لایه خاک برای 
نله کمایک ی کت اما جن وواناق کار ماش اعهت 
حذف اثرات نامطلوب مرز» طول مدل ۱۵-۱۰ برابر عمق 
لایه خاک پيشنهاد می‌شود. چنان چه ابعاد مدل کمتر انتخاب 
شود. مقادیر نیروی محوری بیشتر از مقدار واقعی برآورد 
خواهد شد (بین ۱/۵ تا ۲ برابر) و طراحی بر اساس آن 
ره صرق وه بر خاتان فک انس که فان 
طول مدل کمتر (۶-0 برابر عمق لایه خاک). مشکلی در 


برآورد مقادیر لنگر خمشی ایجاد نخواهد کرد. 


سپاسگزاری 
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۰ 0۰7۱ ,24 ۷۲۵۱ ,۱۱8۵۱6۵۳8۵ ۷۱ وه م۱۱۵۵ ۵ ام ول ۱0۵60 و ولزهگ عقانجهیین طعباهعطا مصتانه۳ 
0۰ 861-879 

۷ عمط ممتاهمه۲۶0۵ ممنتام‌بیک تنج موعبمک گم مصتام۱۷]۵0 لهعزرمصصن اه رطفتامت0 .ه رتطهرصم2 ۷۲۵۲)۵2۵۷1 ,۷۲ [12] 
۰ ,122-132 ۵0۰ ر2 .10 و11 ۷۵۱ ,۳۵۵6۳8۵ ۲) ۵0 ال ۲۱۱۸۵۱۱۵۱۵۸۵ کاصعحصاصده۴ 

تاو ۲ گه مصتاه۵0 لهمزمصاط آمممزمممد ۳۲۳۷۵ بتصتطومطک ,۱۷ مطفتامعم9 هر متصهرصه2 ۱۷۲۵۲۲222۷1 ,1۷۲ [13] 
۰ ,۹8۱66۳۱۱8۵ 0۱۱۱۱۵6 وه م۱۱۵6( ۵۱ ,660۵8۵6 مها ٩۱620۷‏ بصن فصصجنا طاتوظ! طمبامعطا ممتاهعد0۵ظ 
0۰ ,60-71 .00 ,88 

۲۵ ۴5۵۷۵۲9۵ عطا من )حمصصلهطامصط ۵1 مامع]ظ۲ مطا ۵۶ مملاهمناوع۷ ۱ [2ممصصا بتل‌تهصصصه۷ 12 ۷۲۰ مطاویاممو ه رژتمع/ظ ۷۲۰ [14] 
۰ ,113 ,۷۵۱ ,6601۵61076 070 0۳۱۵۱۵۲ 1۵61۵00109 طمتاهصتا۲0 

عصتالناد۲ عوععم؟ ممعباعه ممتامهتعاص1 ۵۶ عصتام۱۷۵0 موتفتتاصمی رما ۷۷۲۰۷۲۰ ,ععنا .2.1 رطم۵012۵0 بخ ر22125ظ ,۱۷۲۰۳۲ [15] 
۰ , .131-164 .00 ,65 ,۷۵1 ب۱۱26۳1۳8ع۱ 1 )0 0۳0 ۱۱۳۸۵/۵6۵ 50۱ ," امصصه 1 220 


سال سی و ششم شماره یک, ۱۶۰۲ نشریه مهندسی عمران فردوسی 


۳ مدلسازی عددی سه بعدی اند رکن شگسل شیب‌لغز و تونل... 


ومصتاهمز۲ ویاممتاومت 0عتتباظ ۵۶ فصمتام۱۷۵0 معبتتممی. پدلمله .و رطمع02722صهلهن0) هر بتل۷]۵۲۵ ,۸۷۲ رتصقطزم۳۴ ,1۷۲ [16] 
۰ ,659-670 00۰ ,49 ,۷۵۱ ,0۵ 6601660171661 ۵۱۵۱۵ بر عصتالناه۲ ۵۷۵۲96 ما 0م6عزطانا 

۵1 1 15مصصت ] «مالقطه لمیمعصوعو گم عصتام۱۵0 ام)ممصصتنهمند۳۴. رطمل0272۵صهاهدان به رتطعقهاطاه .1 رتصدنک ۸۰ [17] 
۰ ,108-119 00۰ ,51 ,۷۵۱ ,766107101088 6۵6۵ ۳۱۵6۲8۵۵۱۵ ۵0۳0۵0 1۱۱۳۱۵۱۱۳۸۵ فالنخ۲ تححصتمآظ بزح ۲1661060 

۵( 0۶ ,۱۷۵008 ]زا۳۲۳2 و2۷۵۲ ۵ ۵0ا0عزهاناو ومصتاهمزظ 0متتباظ 0۶ ولو امصم ردامنز ۸ ررصفطاوه:۲ به بتطوم1 .و [18] 
۰ ,930-940 00۰ ,31 ۷۵۰ ,1۱۱8۱۱62۳1۳8 ۵۲1۵۵6 070 

۵۵5 9601 60اظ ۵۶ ۵2۷10۲ظ ۵۶ ولواقمم اومترممصاظ. رتمک 92112727200 .2 مطفتامتمو به مطم۳]۸۵8۵822۵0 ۷۲۰ [19] 
,1-86 ۱۷۵۷ ب(]0 ماع( ۲ ۲۱۱۵۱۱۵2 از 07۵ و۵۱8۵ ۷۵۱۱۵۱۵۵ (6۸ , فالننج۲ ونله رذن مه 60]640زهانا5 
2014(۰) 

۵ ۵۲۵۷۹0 2110۷ ۵۶ ولو امصخض ممهممیع]۱۱۶ ماتصل اهط10وصمصهز(-۵6ع6:ظ ۳1 رفهزنا .ن1 راتعط22 ,۷۲ رهتصند‌طازصمگ۴ ۱۷۲۰ [20] 
اب اه اه ۵0 ۲۲۵/۵۵۵ 0۵0 عتامجصنحنجظ ظ :۴216 تقصصهماا ه ب۵ الناد۲ ووععم۴ ه عصلوومین) وامصصت [ 
0۰ ,998-1011 .۵0 ,14 ۷۵۱۰ ,۳۲۱۱۵11668 

۰ ,1021102۵1 16۵۲5 2017.01 صم1و۵ ۷ ۲۱2۳5 .عظ1 رفتحه۲۱ [21] 

9۱۱۵ ,60 طمناهعطا ممتاهعدم۱۵۵ظ امین ]۲2 ملفتمطتد۴. بلععو .ظ.۱ پتالن) .وبا بعع8 .ط رهظ .1.1۳ [22] 
4 ,343-961 .00 و3 .20 ,120 ۷۵۱۰ ,۵66108 66016611۱6۵ 

نوی نید وم فعتاتهم۶۵ظ رم لهعنصمطمع۱ 0۴ )۳۲160 رلمه2تتصطاهی نک رطفتامتمو ظ رتصدهزصم2 ۷]۵۲)۵2۵۷1 ,1۷۲ [23] 
2012(۰) راهع ۳0۳ بصمطفنا موه او ون 0۵0۵ ۲۷۷۵۰۵ 7 <وانمگ تمآمهن عمط حمتاععهع۳۸۵ 

عال 0 ۳۱۵2۵۲0 بابک ]انیج۲ 10 ومعاودمع]۱۷ ممتاهع۷۲۱)1 مطلاهع تافه 12۷ بطونامع90 به رتصدزصم2 ۱۷۲۵۲۲۵22۷1 1۷۲۰ ,۳۵211 .۵ [24] 
بطق[ بآ ۵ 0008 ۱۵۱۱۵۵۵۸۵ 3۳ عصناه۱۷۲۵0 امعتتمجصه آظ ره معصناهمزظ 0عتتباظ مه عصتاان۴ منا٩‏ 


(صه۴۵۶ ع) ب(2023) ,17-19 000۵06۲ 


نشریه مهندسی عمران فردوسی سال سی و ششم شماره یک. ۱۶۰۲ 


